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We consider the possibility of quantum phase transitions in the ground state of triplet supercon- 
ductors where particle density is the tunning parameter. For definiteness, we focus on the case of 
one band quasi-one-dimensional triplet superconductors but many of our conclusions regarding the 
nature of the transition are quite general. Within the functional integral formulation, we calculate 
the electronic compressibility and superfiuid density tensor as a function of the particle density for 
various triplet order parameter symmetries and find that these quantities are non-analytic when a 
critical value of the particle density is reached. 

PACS numbers: 74.70.Kn,73.43.Nq 



Triplet superconductivity is a very rare, but very rich 
phenomenon in condensed matter physics. The very few 
confirmed examples in nature include strontium ruthen- 
ate [l| (a Ruthenium oxide) and the Bechgaard salt 
(TMTSF) 2 PF 6 2] (an organic molecular compound). 
Because the full confirmation of triplet superconductiv- 
ity in solids has occurred only over the last few years, 
these lattice systems are not yet as fully studied theo- 
retically and experimentally as 3 He, their neutral liquid 
superfiuid counterpart Q . Unlike 3 He, these systems are 
lattice charged superfluids, and their order parameters 
are intimatelly related to the lattice periodicity as in d- 
wave high critical temperature (T c ) superconductors p|. 

It is known experimentally that electronic properties 
and ground states of cuprate superconductors (d-wave 
singlet) 0| and Strontium Ruthenate (p-wave triplet) || 
are very sensitive to chemical doping while these prop- 
erties for (TMTSF) 2 PF 6 (p-wave triplet) are very 
sensitive to both external pressure and chemical dop- 
ing. However, recent experiments have demonstrated 
that it is possible to change the carrier density electro- 
statically in cuprate superconductors 6] and amorphous 
Bismuth without the introduction of additional dis- 
order that often occurs for chemical doping. It is likely 
that similar electrostatic techniques will be developed 
for use in other superconductors like Strontium Ruthen- 
ate or (TMTSF) 2 X (with X = C10 4 , PF 6 ). Since it may 
be possible in the near future to tune the carrier con- 
centration with the use of field effect techniques, some 
important theoretical questions concerning these three 
types of systems may soon receive an experimental an- 
swer. For instance, are there quantum critical points sep- 
arating magnetic and superconducting order as a func- 
tion of particle density? Or, are there quantum critical 
points within the superconducting phase as a function 
of particle density? In Fig. ^ we show the zero tempera- 
ture density versus interaction phase diagram indicating 
the existence of quantum critical lines, where the order 
paramater does not change symmetry but the ground 



state topology changes. The elementary excitation spec- 
trum also changes from gapless to gapped, and at finite 
temperatures there are three distinct regions the Fermi 
liquid, the pseudo-gap, and the Bose liquid regions. 

In anticipation of experimental efforts, we propose 
to study the possible existence of topological quantum 
phase transitions in lattice triplet superconductors as 
a function of particle density. We focus on the spe- 
cific case of the (TMTSF) 2 X family, because they are 
a single-band triplet superconductors, unlike Strontium 
Ruthenate where three bands may be necessary to de- 
scribe triplet superconductivity 0. For this purpose, we 
study single band quasi-one-dimensional systems in an 
orthorhombic lattice with dispersion 

ek = — t x cos{k x a) — t y cos(k y b) — t z cos(fc z c), (1) 

where t x ^> t y ^> t z . We work with the Hamiltonian 
H = Hkin+Hinti where the kinetic energy part is Hu n — 
J2k, a ^kV , k jQ V'k,Q, with £ k = e k -/i, where fi may include 
the Hartree shift. The interaction part is 

H ™t = \ J2 £ K l /3 75 (k,k')6i /3 (k,q)6 7 5(k',q) (2) 



kk'q afi-yS 



with &^(k,q) = 



t 



k+q/2,Q^k+q/2,/5 



, where the labels a, 



f3, 7 and 5 are spin indices and the labels k, k' and q 
represent linear momenta. We use units where h~ks = 
1. In the case of weak spin-orbit coupling and triplet 
pairing, the model interaction tensor can be chosen to 
be 

7 a/37 *(k,k') = -|VtMk,k')<Mk)<^(k')r Q/375 , (3) 



where T 



vL/2 with v Q/3 = {iaa y ) a p. V r 



is a prefactor with dimensions of energy which char- 
acterizes a given symmetry. Furthermore, the term 
/ir(k, k')0r(k)<^p(k') contains the momentum and sym- 
metry dependence of the interaction of the irreducible 
representation T with basis function 0r(k) and 0p(k') 
representative of the orthorhombic group (-D 2 /i). 
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We write down the partition function [8| 
J V lip', t/j] exp [— S] with the action 

For 



(a) 



z 
s 

simplicity, we take /Vr(k, k') = 1. We introduce the 
vector order parameter for triplet superconductivity 
through the Hubbard-Stratanovich transformation and 
integrate out the fermions to obtain the effective action 



S e ff = Qg 



dr Ey-^ lTrln 



(3M 



(4) 



where Q G = J ^r £ q 4 Pj(q, T)2Mq, r)/|F r | and M is 
the matrix 



M = 



[d r + 6k] 5 klk2 5 a/3 A(ki,k 2 ,r) 
At(ki,k 2 ,r) [d T - 6J S^^ap 



where A = - £\ r (^^)'D i (k 1 + k 2 , t)v q/3i1 . At the 
saddle point, £>j(ki-|-k 2 , r) is taken to be r independent, 
and to have total Fermion center of mass momentum 
ki + k 2 = 0. Thus, 2>i(ki + k 2 ,r) = ^Ar^+ks.o + 
*X\(ki+k a ,T). 

The order parameter equation is obtained from the 
stationary condition 6Sj£) /5T>j — 0, where S^S is the 
saddle point action 



1 



E 

k 



y r ||0r(k)| 2 tanh(/3£; k /2)/2£; k 



(5) 



where E^ = y6k + Ar| 2 |0r(k)| 2 , corresponds to the 
quasiparticle excitation energy. The number equation is 
obtained from N = —dil/dfi, where /3f2 = — InZ is the 
thermodynamic potential, 



N = N Q +N Ruct , 



(6) 



where N = X) k n k, and rik = [1 — £ktanh(/3i?k/2)/-Ek] 
is the momentum distribution. The additional term 
A'fluct = — Sf2fl uc t/9A<, where r2fl uc t are Gaussian fluctu- 
ations to saddle point f2n. These two equations must be 
solved self-consistently, and quite generally they are cor- 
rect even in the strong coupling (or low density) regime 
provided that T < T c 8], where 2V fluct ~ T 4 for all 
couplings |9|. 

The saddle point vector T>^ = rji Ar is related 
to the standard d- vector via the relation <ij(k) = 

Ek' Pf'frW' In the D 2 h point group all represen- 
tations are one dimensional and non-degenerate flfj| . 
which means that the d-vector in momentum space for 
unitary triplet states in the weak spin-orbit coupling 
limit is characterized by one of the four states: (1) 
3 A lu (a), with d(k) = fjA f ^XYZ («f xyz » state); (2) 
3 B lu {a) 1 with d(k) = r,A p ^Z ("p z » state); (3) 3 B 2u (a), 
with d(k) = f)A Py Y ("p y n state); (4) 3 B 3u (a), with 
d(k) = -qAp^X ("p x " state). Since, the Fermi surface 
can touch the Brillouin zone boundaries the functions 
X, Y, and Z need to be periodic and can be chosen 
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FIG. 1: Phase diagrams for (a) f xyz (b) p x symmetries based 
on Fermi surface connectivity and excitation spectrum. The 
phase diagrams are symmetric (not shown) around half-filling 
(N = 1) due to particle-hole symmetry. The small cross in- 
dicates the parameters compatible to the known triplet su- 
perconductor (TMTSF) 2 PF 6 N = 0.5 and V p Jt x = 0.5531 
for p x , N = 0.5 and Vf xy Jt x = 2.1685. 



to be X = sm(k x a), Y = sin(k y b), and Z = sin(fc z c). 
The unit vector i) defines the direction of d(k). From 
here on we scale all energies by t x . The parameters used 
are t x = 1.0000, t y = 0.2144 and t z = 0.0083. It will 
be easier to discuss the properties of the Fermi surface 
and quasiparticle excitation spectrum by considering the 
chemical potential (x instead of the filling factor N. Con- 
sider the "normal state" Fermi surface defined in the 
first Brillouin zone (BZ) by 6c = 0, keeping in mind the 
periodicity in fc-space. Let us define the following char- 
acteristic values p,\ = t x + t y + t z , p\ = t x + t y — t z , 
/ig = t x —t y +t z , p\ = t x — t y — t z . For p < p* the chem- 
ical potential is below the bottom of the band and there 
is no Fermi surface. For pi < p < p\ the Fermi surface 
consists of one connected sheet contained entirely in the 
first BZ and is topologically equivalent to a sphere of 
genus zero. For p^ < fi < /I3, the Fermi surface is one 
connected sheet touching the edges of the BZ on the 
planes k z = ±n equivalent to a torus of genus one. For 
p% < fi < pi, the Fermi surface is one connected sheet 
touching the edges of the BZ on the planes k z = ±n 
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and k y = ±7r equivalent to a torus of genus two. Lastly 
for \i\ < /j,, the Fermi surface consists of two discon- 
nected sheets each of which touch the BZ boundary and 
is equivalent to two disconnected tori each of genus one. 

For the superconducting state, the intersection of the 
Fermi surface and order parameter nodes constitute the 
loci of gapless quasiparticle excitations. For the Pi sym- 
metry (where i is x, y, or z) the order parameter nodes 
are on the planes fcj = 0, ±7r. For the f xyz symmetry, 
the order parameter nodes are on the union of the planes 
fej = 0, ±7r, where i is x, y, and z. For all symmetries, 
the quasiparticle excitations are fully gapped for li < p\ 
since there is no Fermi surface. For ll\ < fi < n\, the 
order parameter nodes for all symmetries intersect the 
Fermi surface and hence quasiparticle excitations are 
gapless. For /x| < /i the Fermi surface splits into two 
sheets that separate along k x so that the p x nodes no 
longer intersect the Fermi surface opening a gap in the 
quasiparticle excitation spectrum. However, the f xyz , 
p z , and p y nodes still intersect the Fermi surface so quasi- 
particle excitations remain gapless. 

We plot representative phase diagrams for the f xyz 
and p x symmetries in Fig. ^ From the discussion above, 
the p y and p z phase diagrams are qualitatively similar 
to the f xyz phase diagrams. There are three distinct 
phases characterized by Fermi surface connectivity and 
quasiparticle excitation spectrum: (1) no Fermi surface 
and fully gapped E(k) for all symmetries (fi < /Lt*), (2) 
one sheet Fermi surface and gapless for all symmetries 
(fj,* < fi < (3a) two sheet Fermi surface and gapless 
for fxyz, Pz, Py (fJ't < p) (3b) two sheet Fermi surface and 
fully gapped for p x (/x| < //). In addition, the (2) phase 
splits into three regions using the finer classification of 
Fermi surface topological genus: (2i) genus zero (/x* < 
/it < (2ii) genus one (fx^ < P < fJ%), (2iii) genus two 
(/I3 < p < There are several qualitative features of 
interest in the phase diagrams. For a fixed density and 
all symmetries, the chemical potential does not go below 
the bottom of the band (p < /j,*) until a critical coupling 
is reached. This critical coupling is Vf x z /t x = 14.8021, 
V p Jt x = 1.8150, V Py /t x = 2.3952, V p Jt x = 3.5052 for 
N = 0. In addition, consider the approach to the strict 
one-dimensional limit (t y ,t z — > 0). If t z — > 0, the ii\ 
(/ig) boundary merges with li\ {p%) leaving only the (2ii) 
region. If in addition t y — > 0, the li\ boundary merges 
with ii\ producing only one boundary between the (1) 
and (3) phases. 

We now turn our attention to thermodynamic quan- 
tities that provide signatures of the topological changes 
discussed above. In the following calculations we fix the 
interaction strength to be Vf x z /t x — 91.7202, V Pz /t x — 
11.8747, V p Jt x = 11.7781, V Vvs /t x = 10.8297, 'which 
forces Li = fi* at N = 0.5. In contrast, fi = ll% at 
N = 0.706, N = 0.674, N = 0.672, N = 0.719 for the 
fxyz, Pz, Py, Px symmetries, respectively. The T = 




FIG. 2: Filling factor (TV) dependence of the dimensionless 
electronic compressibility n\t x \. The inset shows the region 
JV = 0.4- 0.8. 



electronic compressibility n = N 2 (dN / dp)x,v is 



2E k \ El 



(7) 



which we plot in Fig. [2] Although the compressbil- 
ity does not formally diverge, there are clear anoma- 
lies (non-analyticities) when the Fermi surface topol- 
ogy and quasiparticle excitation spectrum change. As 
N increases, y, increases and crosses the boundaries [i* 
where the first derivative of k decreases discontinuously. 
Within the p-symmetries, the magnitude of this jump is 
largest for p x . The f xyz symmetry, with the presence 
of double or even triple nodes compared to the single 
nodes for the p symmetries, has the largest jumps in the 
derivative of n clearly identified as the four cusps in Fig. 
12 These non-analyticities in k at T — are indicative 
of a quantum phase transition. At finite temperatures 
the cusps in k are smeared-out, but clear peaks are still 
present so long as one remains in the quantum critical re- 
gion. The measurement of the electronic compressibility 
may be achieved in a field effect geometry 0, @ through 
the relation 



n = VC d /Q 2 , 



(8) 



where Ca = [dQ /dV e ]T,v is the differential capacitance, 
V e is the applied voltage, Q is the absolute value of the 
total charge of carriers, and V is the sample volume. 

Next, we analyse the effect of phase fluctuations given 
by the action 

AS = 1 12 tt^K) 2 + Pum] ( 9 ) 

where A — N 2 k/V is proportional to k and 



Pij = y 12 [«k9»9j^k] , 



(10) 



is the superfluid density. Here, is the momentum 
distribution. In Fig. we show the N dependence of 
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FIG. 3: Filling factor (N) dependence of the zz component 
of the dimensionless superfluid density tensor Vop zz /a 2 z \t x | for 
various symmetries. Vo (a z ) is the unit cell volume (length). 

the zz component of pij. In the case of the Dih group 
only diagonal components pa exist, but they are highly 
anisotropic due to the quasi-one-dimensionality of £k- 

From Eq. I1UI it is clear that the zero-temperature 
superfluid density is the curvature of the dispersion 
weighted by the momentum distribution. We find the xx 
component is a monotonically increasing function of N 
which is best understood as a consequence of quasi-one- 
dimensionality. Neglecting the curvature of the Fermi 
surface due to finite t x , t y , as N increases, the Fermi 
surface encloses a monotonically increasing pocket of k- 
space around k x = 0. In addition, below half filling 
the curvature of is of the same sign in this regime so 
there is little cancellation between different regions of k- 
space. We find p yy is also smooth but peaks in the region 
bounded by N « 0.5 — 0.8 corresponding to fi « /x^ — p,\. 
No longer neglecting the curvature of the Fermi surface 
due to finite t y , we see that for N « 0.5 — 0.8, the Fermi 
surface encloses an increasing pocket of fc-space around 
k y = up to the edge of the BZ. Since the curvature 
of £fc changes sign at k y — tt/2, cancellations between 
different regions of fc-space eventually occur. The zz 
component is the most interesting as it exhibits clear 
anomalies such as those seen in k. The Fermi surface 
varies rapidly along the k z direction as a function of N 
so that no simple analysis in terms of contributing re- 
gions of fc-space is possible. It appears that the p zz is a 
more direct probe of the anomalies seen in k with clear 
kinks as a function of N. These non-analyticities again 
indicate the existence of a quantum phase transition as 
/i is tuned below the bottom of the band. 

Lastly, it is easy to extract from AS the phase-only 
collective mode frequencies via the substitution ioj n — ► 
uj+id. In this case w(q) = yj (c x q% + c y q y + c z q z ), where 

cl = Pxx/A c 2 y = p yy / 'A, and c\ = p zz /A. These col- 
lective mode frequencies, also show similar anomalies 
to those of k and p^j, and can be used to character- 
ize the quantum phase transition as well [9j . Notice that 



w(q) is anisotropic reflecting the orthorhombic structure. 
For instance, u)(q x ,0, 0) = c x q x , Lo(0,q y ,0) = c y q y , and 
w(0, 0, q z ) = c z q z , where Cj is the speed of sound along 
the i-th direction. However, these modes may be plas- 
monized in a charged superfluid. 

In summary, we studied possible quantum phase tran- 
sitions in triplet superconductors, as the density of car- 
riers is changed, provided that the Cooper pairing in- 
teraction is sufficiently attractive. For organic quasi- 
one-dimensional conductors (Bechgaard salts) only one 
quantum phase transition may be accessible experimen- 
tally as the interaction strength is too weak to cross 
both phase boundaries as a functions of filling factor (See 
Fig.nj. However, these quantum phase transitions may 
be possible in optical lattices where the attractive in- 
teraction may be changed via Feshbach resonances, thus 
allowing the system to cross the phase boundary sepa- 
rating the weak and strong interaction regimes [llj. In 
order to identify these quantum phase transitions (QPT) 
where the symmetry of the order parameter does not 
change, we classified Fermi surface topologies and exci- 
tation spectrum properties as a function of filling fac- 
tor for weak spin-orbit coupling symmetries f xyz , p x , 
p y and p z . We then related non-analyticities in the elec- 
tronic compressibility and superfluid density to quantum 
phase transitions between various phases. We would like 
to thank NSF (Grant No. DMR-0304380) and NDSEG 
for support. 
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